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We describe a class of the singular solutions to the multicomponent analogs of the Lame 
equation, arising as equations of motion of the elliptic Calogero-Moser systems of particles 
carrying spin 1/2. At special value of the coupling constant we propose the ansatz which 
allows one to get mcromorphic solutions with two arbitrary parameters. They are quantized 
upon the requirement of the regularity of the wave function on the hyperplanes at which 
particles meet and imposing periodic boundary conditions. We find also the extra integrals of 
motion for three-particle systems which commute with the Hamiltonian for arbitrary values 
of the coupling constant. 
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This letter is devoted to the problem of finding solutions to the matrix equation which arises 
naturally in the theory of quantum Calogero-Moser A^-particle systems [l|, y| , 
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ij = Eip, 



(1) 



where {xj} are coordinates of the particles, Xjk = Xj — Xk, {sj} are their spins, ip depends on 
{xj} and {sj}, a € R+, {Sjk} permute Sj and Sk and form the SU{n) spin representation of the 
permutation group Sn^ and p{x) is the Weierstrass elliptic function with two arbitrary periods 
2uji, 2u;2 with only restriction Im(u;2'-^r ) t^ 0- ^^ shah suppose for definiteness that uji is real. 
At A^ = 2, {Sjk} = 1, one arrives at the usual Lame equation, but the matrix problem is highly 
nontrivial even in the case of A^ ^ 3, n = 2 (particles carry spin 1/2) which we shall mainly discuss. 
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Various limits of the equation ([T]) were extensively studied in a lot of papers 
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la . llSl . lid ] and references therein. In [3], Olshanetsky and Perelomov found the 
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connection between ([T|) and the roots of ^at-i algebra, and proved the integrabihty of quantum 
spinless Calogero-Moser systems (i.e. {Sjk} = 1) but did not find any way of determining ip. The 
first results for spinless particles were obtained in [41] where the explicit form of ip was found for 
o = 1, A^ = 3, and algebraic structure of the manifold containing all -(/''s was described for all 
a € N, and all N. Later on, overcomplicated meromorphic solution for all a E N, A^ was found 
by Felder and Varchenko [5]. In [Q], the author obtained explicit formulas for spinless case for all 
a € R+ in the form of infinite series. 

so intensively studied [71, |8|, M, [l 



The trigonometric limit of ([T]) (|w2| -^ oo) has been a 
121 . I13I . Ilq for the case of particles carrying spin. In [T], la. 
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111 ], the solutions were found as 



the spin generalization of the Jack polynomials which provide adequate description in this limit 
of the eigenvalues problem for spinless case and arbitrary a E R+. The corresponding symmetry 
responsible for these exact results was found to be the y[SU{n)^ Yangian algebra 12, llJ, 114] • 
The authors of the paper [16|] have considered elliptic spin case but they have found the effect of 
(quasi)exact solvability only for the case of the BCn root systems, and ([1]) does not fall into their 
classes. The quantum Lax pair, i.e. the solution to the equation [H,L] = [L,M], L, M being 
(A^ X A^)-matrices, was also mentioned in [15] for the elliptic case, but it was constructed without 
any dependence on a spectral parameter. Moreover, in the elliptic case even the existence of the 
Lax relation [H, L] = [L, M] does not give the integrals of motion in the form of Ij = Y2j^ i{U)ki 
since the M matrix does not obey the "sum-to-zero" conditions X]i=i ^jk = Sfc=i ^jk = 0. 

In this situation, to our mind, every analytical results to the solutions of the elliptic matrix 
Schrodinger equation ([T|) are of value, even for some restrictions for a, A^ and n. In what follows, 
we put A^ = 3 (three-particle case) and consider at first the question of the integrability of the 
problem defined by ([1]) for arbitrary a and n. For the spinless case, it is known [3] that there is 
the operator 
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(2) 



commuting with H in this case, where a is an arbitrary "spectral" parameter which arises in the 
Lax-matrix approach. One can try to generalize this structure for the spin case as 



a.-^^^^^^+^ !:«(«+ ^.'^)(P(-.'=)-p(«))^' (3) 
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but direct calculation shows that (l3|) does not commute with H and the term of higher order in 
permutations must be added. To gain some intuition for constructing the proper operator, we 



consider the "freezing trick" (a -^ cxd) by which the integrals of motion for quantum elhptic spin 
chain [l7| should be obtained, and try to add the analogs of these integrals to the ansatz ([3]), 

^^^^" = dxjx2dxs + i E «(« + S,k){pix,,) - p{a)) — 

3 

+ A ^ f{xjk)f{xki)f{xij)SjkSki, (4) 

where the function f{xjk) is analogous to the elements of the Lax matrix for quantum spin chain 
situation, 

/(.,,) = !lpi±^ (5) 

a{xjk)cr{a) 

Here a{x) is the Weierstrass sigma function, ^^ log cr(x) = — p(x), a{x) ~ x + O(x^) as x — > 0, 
and A is some parameter which should be determined by the commutativity condition, 

[-H , JJspinl ^ U. 

Our computations shows that the ansatz (jl]) is indeed correct and 

A = aV3. (6) 

Moreover, the function 

^jki = f{xjk)f{xki)f{xij) (7) 

is simplified drastically when being considered as the elliptic function of the spectral parameter 
a, having a pole of third order as a — > (and no other singularities on the torus C/F, F = 
2Za;i + 2Zu;2)- It can be written as 

^jki = - 2P'(") + '^jkiP{a) + ifjkh (8) 

where the functions tpjki and (pjki do not depend on a and can be found by calculating the coeffi- 
cients of the Laurent series for ^ ai a —> Q, 

^jki = C{xjk) + C{xki) + Cixij), (9) 

^jki = --^{p'i^jk) + p'ixki) + p'ixij) + 2[C{xjk) + Cixki) + C(a;/j)]^|- (10) 

The formula (jlOp is obtained with the use of the well-known relation 

p{Xjk) + P{Xkl) + p{xij) = [C{Xjk) + Cixki) + C{xij)] , 



where ((x) is the usual Weierstrass zeta function, C'i^) = ~p{x)- Hence we found two independent 
integrals of motion from (H HlOp due to the arbitrariness of a spectral parameter a, 
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j<k<l ^ '' j^k^l^j j^k^l^j 
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^2=2 E ^^^^~Z E ^JklSjkSkl- (12) 

The formula for ^2 is especially simple: it resembles the total momentum (and coincides with it 
as {Sjk} = !)• In the trigonometric limit, it can be expressed through the scalar product of the 
Yangian generator and total spin which (in this limit only!) both commute with H. We confirmed 
also, by direct computation of [j7i,^2]) that these operators mutually commute and form with H 
and total momentum the commutative ring for all values of the parameter a. 

Let us now construct the explicit solutions of ([T]) for the simplest nontrivial case of three 
particles carrying spin 1/2. When all spins aligned up or down, we have the situation analogous to 
the spinless case [J]. The nontrivial form of the wave function ip arising for the states with total 
spin S = 1/2 is as follows 

V^(xi,a;2,X3)=^(x)|Tn> + i?(x)inT) + C(x)UTT), A + B + C = 0. (13) 

The operators {a{a + Sjk)} act on the spin pairs in the states (IT), (Ti + it) as a{a + 1) and for 
(Ti ~ it) as a{a— 1). If a is chosen as positive integer, there are singularities in the spinless case in 
the form of poles, {xj — Xk)~"' as Xjk — > 0. It is natural to expect that in the case of particles with 
spin at least some solutions to ([T|) have the similar structure, i.e. A, B, C have singular behavior 
as Xjk — > in the form of poles. The equation ([1]) reads in the component form after substituting 
([ED as 

(l ^ Q2 \ 3 

oEtT^+^M-E P(^jfc)^ - [P(a^l2)^ + P{XZI)C + P(X23)5] = 0, (14) 
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oEtT^ + ^^-E ^(^J-fc)^ - [pi^i2)C + p{xu)B + p{x2i)A] = 0, 
^ 3=1 "^""i / ]>k 

I '^ Q'i \ 3 

2^^^^r~Y^ P(^Jfc)C - [p{xi2)B + p{x^i)A + p(x23)C] = 0. 



(15) 
(16) 



i=i ■? / j>k 

Let us introduce the notation 



y(x) = A(x) - S(x), Z(x) = A(x) - C7(x) (17) 



and deduct p5|) and (|T6|) from (|M|) . Under the condition p3j) . it is easy to see that the system 
(fT3HT6]) is equivalent to two coupled equations for Y{x) and Z{x), 

2 ^ a^ + ^ r " ^^^12^^^ + ^) " ^^^ " ^)P(^31) = 0, (18) 

\ j=l j / 

(^E^ + ^)^-^(^i2)(^ + ^)-(2^->^)p(x23) = 0. (19) 

Since the "potentials" here are double periodic, it might be expected that the solutions to 
(jl8|19p are quasiperiodic, acquiring the same Bloch factors under the shifts of the arguments by 
the periods 2u;i, 2lo2 of the Weierstrass functions. According to (|18p . l^(x) has a simple pole at 
3^31 -^ 0, the same is for ^(x) as X23 -^ 0. The analysis of limits X23 -^ for (J18p and X31 — > for 
(J19p shows that the left-hand sides of (|18p and (J19p are regular at these conditions. And finally, if 
y(x) ^ ^(x) as j;i2 ^ 0, there should be a simple pole singularity of these functions in this limit. 
Combining all these properties, we come to the ansatz for Y and Z in the form 

^, . ,o-(/xi2)ct(xi2 + Ai2)o-(2;3i +A31) , t. N /onA 

y(x) = h exp(fcixi + ^2X2 + ftsxs), (20) 

0-(j;i2)cr(X3i) 
ry, . ,0-(Al2)0-(2;i2+^12)0-(x23 + /i23) ,, , , , , ^ ln^\ 

Z{x) = b exp(/cia;i + ^2X2 + A;3a;3), (21) 

a{xi2)cr{x23) 

where a{x) is the sigma Weierstrass function defined above, 6, {A;j},Ai2, A31, /ii2, fJ-23 are some 
parameters. The Bloch factors for (j20p and ()2ip are equal if and only if 



/^12 — A12 — A31, /i23 — — A31. (22) 

These expressions look rather asymmetric in {A,B,C}, but the symmetry becomes evident with 
the use of the remarkable identity 



y(x) - Z(x) = C(x) - B(x) 

cr(A3i)cr(2;23 - Ai2)cr(x3i - A12 + A31) 



exp(A;ixi + A;2X2 + /i:32;3), (23) 



cr{x23)cr{x3i) 

which is valid for all values of the parameters A12 and A31 and coordinates {xjk}- Some long but not 
too tedious calculations show that (i20]) and (f2T]l under the condition ([22]) indeed give the solutions 
to the system (|18|19p if the following restrictions to the parameters {A}, {k} take place, 

ki-k2 = C(A3i - A12) - C(Ai2), (24) 

A;2-A;3 = C(A3i) + C(Ai2), (25) 



{ki,k2,h) 


Al2 


A31 


-0(x,s,k,Ai2,A3i) 


{ki,k2,k3) 


Al2 


A31 


-0(x,s,k,Ai2,A3i) 


{k2,ki,k3) 


— (Ai2 — A31) 


A31 


-ni2V'(x,s,k,Ai2,A3i) 


{ki,k3,k2) 


-A31 


— A12 


-n23'0(x,S, k, Ai2,A3i) 


{k3,k2,ki) 


A12 


A12 — A31 


-n3iV'(x,s,k,Ai2,A3i) 


{k2,k3,ki) 


-A31 


A12 — A31 


ni2n23V'(x,S,k,Ai2,A3i) 


{k3,ki,k2) 


— (A12 — A31) 


— A12 


ni3n32-0(x,S,k,Ai2,A3i) 



TABLE I: Relation of the transformations of the parameters {k} and {A} of the eigenfunction tp that obey 
and the action of permutation operators of particles {Hij} on ip. 



where Ci^) is the zeta Weierstrass function defined above. To get (|24|25p . we used the formula 

>./ N ^/ N ^/ N >/ X crlx + y)a(y + z)a(z + x) 

C{x) + C y + Ciz) -ax + y + z)= . '{ )\ /2 ^ V 

a{x)a{y)a{z)a{x + y + z) 

As for the corresponding eigenvalue, it can be written in very symmetric form, 



E 



--(/Cl + k2 + kif - -(p(Ai2) + p(A3l) + p(Ai2 - A31)), 



(26) 



with two still unspecified parameters A12, A31. Note that (|2UI21I26|) look very similar to the solution 
of the usual Lame equation, 

d?il){x) 



dx"^ 



+ a{a + l)Tp{x) = Eipix), 



in the case of a = 1, where ip{x) ~ exp ( — xC,{a))a{x + a) [cr{x)] (the Hermite (1872) solution), 
and E = — p(a). We also found by some long calculation that (|2UI2ip form also the eigenfunctions 
of the operators J7i, J2 (|llll2p at a = 1, with the eigenvalues 

Jl = ^{k,+k2 + hf - ^^^t'^^' (C(Al2) - C(A3l) - C(Al2 - A31)) 
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14(C(Ai2) - C(A3i) - C(Ai2 - A31))' + 9(p'(Ai2) - p'(A3i) - p'(Ai2 - A31)) , (27) 



j2 = C(Ai2)-C(A3i)-C(Ai2-A3i), 



(28) 



also with arbitrary values of the parameters {A}. 

Since the relation (|26|) is invariant under all the permutations of the indices of {k}, the complete 
symmetrization of (J13p on permutations of particles also gives an eigenfunction of H, we named 
it '0o(x, s, k, A12, A31), if {k} and {A} obey (|24H26p . The complete set of transformations of the 
parameters and its relation to permutation operators of particles is contained in Table [H This 



eigenfunction ^o(x, s, k, A12, A31) is regular as Xj — rr^ ^ 0. Furthermore, Eq. ([26]) is invariant 
under global sign reversals of {k} and {A} so we have another eigenfunction il^i{x, s, k, A12, A31) also 
symmetric under permutations of particles and regular as Xj — Xk ^ and linearly independent of 
^o(x, s, k, A12, A31), that verifies 

^i(x,s,k, Ai2,A3i) = ^pQ{yL,s, -k, -X12, -X31) = -'0o(-x, s,k, A12, A31). (29) 

The discrete spectrum of the corresponding three-particle system on the real circle 3^1,2,3 G 
Rmod2u;i can be obtained by imposing the periodic boundary conditions 

{ki - k2)uji - C(wi)(Ai2 - 2A31) = ivrZi 

/i,/2GZ. (30) 

{h - h)^i - C('^i)(Ai2 + A31) = ivr/a 

There is a degeneration in the energy, but ipQ and ipi have distinct and opposite eigenvalues 
through the action of J'l and J2 ■ Due to the relation satisfied by the total spin 

1113 11 

-(8)-(8)- = -©-©-, 
222 222' 

there exists an intrinsic degeneration of functions with total spin 1/2. That is the reason why we 
have such a degeneration on the energy level (j26p . 

To conclude, we obtained for the first time the extra integrals of motion for the elliptic Calogero- 
Moser system of three particles with spin (jllll2|) . Thus we proved its complet integrability. It can 
be shown by direct computation with the use of the Liouville theorem that, replacing 3 by 4 in 
(|ll|12p . J'l and J2 are also two mutually commuting integrals of motion for the A^ = 4 case. We 
found two non-trivial meromorphic eigenfunctions, depending on two parameters, for the spin 1/2 
case and coupling constant a = 1 ()20l2ip . 

We can conjecture that the regular eigenfunctions of ([T|) will be totally symmetric functions 
under permutation of particles for all values of N and a as it has been shown for some of its limits 
[13] and for other unrelated spin dynamical models J18l |. 
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